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A procedure is described for the precise nonrelativistic 
evaluation of the dipole polarizabilities of H2 + and D2 + that 
avoids any approximation based on the size of the electron 
mass relative to the nucleus mass. The procedure is con- 
structed so that sum rules may be used to assess the accuracy 
of the calculation. The resulting polarizabilities are consistent 
with experiment within the error bars of the measurements 
and are far more precise than values obtained by other theo- 
retical methods. 

PACS numbers: 33.15.Kr, 33.15.-e, 31.15.Ar 



The separation of nuclear and electronic motion is the 
underlying principle of the theory of molecular struc- 
ture. The theory is challenged by recent measurements 
of Jacobson et al. |j| of the electric dipole polarizabil- 
ities of H2 + and D2 + which have a precision beyond 
that obtained in the Born-Oppenheimer approximation. 
The measurements stimulated the introduction of meth- 
ods l^-H that take into account the diabatic coupling 
omitted in the earlier calculations and they led to polar- 
izabilities that agree with the measured values within the 
combined experimental and theoretical uncertainties. We 
present here new theoretical predictions of much greater 
accuracy which in turn pose a significant challenge to ex- 
periment. The accuracy of our method can be assessed 
by the use of sum rules and we predict nonrelativistically 
the polarizabilities of H2 4 " and D2 + to a precision well 
beyond that achieved by the experiments. The method 
is general and it should be possible to apply it to many- 
electron diatomic molecules. 

Separating out the center of mass motion we may write 
for the Hamiltonian of H2" 1 " or D2" 1 " in an electric field 
F = Fn lying along the Z-axis of the space-fixed frame 



H 



517 V fi I^ + Tm)^ 



V(r,R) + (1 + e)Fh-v, 
(1) 

where R is the vector joining the nuclei, r is the position 
vector of the electron measured from the midpoint of 
R, M is the mass of the proton or deuteron, V(r, R) 
is the electrostatic interaction potential and (1 + e) = 
[1 + (1 + 2M)~ 1 ]. We use atomic units throughout. The 
change in energy of the system for small values of the 
applied field is given by AE = — ^ctdF 2 , where ad is the 
polarizability. Thus if ^^(r, R) is the eigenfunction of 
the unperturbed system with Hamiltonian Hq and Eq is 
the eigenvalue, the polarizability can be written 



where 

{Ho - #o)* (1) (r, R) + (1 + e)n • r¥<°> (r, R) = 0. (3) 

Alternatively ll^ 1 ) can be determined from the stationary 
value of the functional 

J = (* (1) |#o - Bb|* (1) > + 2(1 + e)(* (1) |n • r|*(°>). (4) 

If we write (r, R) as an expansion over some chosen 
basis set ip n (r, R), 



N 



*«(r,R) = ^Q^„(r,R), 



(5) 



n=l 



assumed to diagonalize the unperturbed Hamiltonian Hq 
so that (ipn\Ho\tp n >) = E n S nn >, the polarizability may be 
written 



ai = 2(1+e)2 £ffi!^^!. 

„ E n — E 



(6) 



This expression for the polarizability is stationary with 
respect to first order errors in and is bounded from 
below. 

The completeness of the set ^ n (r, R) can be assessed 
by inspecting other sum rules. Introduce the oscillator 
strength 



/„ = 2[(E n E )/(l + 3k)]|<¥<°>|n ■ r|^)| 2 
and define the sum 



S(p) = J2[(E n -E )/(l 



2M ' 



r 

J n 



71=0 



so that 

a d = (l + e) 2 (l + ^)- 1 -5(-2). 
Then provided the tp n form a complete set, 

s(-i) = §<*( o yi^ )) 

and 

S(0) = 1. 



(7) 
(8) 

(9) 
(10) 

(11) 



a d --2(*( 1 )|(l + e)n-r|*(°)), 



(2) 



The eigenf unctions $'°)(r, R) and )/> n (r,R) can be 
written as sums of products of nuclear and electronic 
wave functions of the form 



1 
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4tt 
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^ A (e,$,O)0 sA (r,i?) XsA (i?), 

(12) 



where (O, <&) are angles specifying the orientation of the 
internuclear axis in the space-fixed frame, N is the total 
angular momentum quantum number, M is the projec- 
tion on to the space-fixed Z-axis, A is the projection of 
the electronic angular momentum on to the internuclear 
axis and D is the rotation matrix |g|. For the ground 
state of H2 + or D2 + , N = M = A = and the electronic 
wave function has symmetry. The perturbed state is 
a superposition of states with N = 1, M — and A = 
and ±1, the electronic wave functions having and Tl u 
symmetry. 

To calculate the matrix elements of the Hamiltonian 
and of the electric dipole operator we transform from 
the space-fixed frame to the body-fixed frame following 
standard procedures [^]]7]]. The nuclear kinetic energy 
operator may be written 



V 2 

v R 



2M 

where H rot is given by 
1 
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(13) 



2MR 2 
1 

2MR 2 



(N - L)^ 

(N 2 + L 2 - N~L+ - N+L- -2A 2 ), (14) 



in which L is the electronic angular momentum and ± 
indicates angular momentum raising and lowering opera- 
tors. These are the operators that couple £ and II states. 
The electronic wave functions for H2 4 " and D2 4 " are sepa- 
rable in prolate spheroidal coordinates and we expressed 
the electronic basis functions <^ q a(i", R) in terms of these. 
The corresponding formulas for the matrix elements of 
H mt are given by Moss and Sadler J§| . A detailed descrip- 
tion of the representation of the nuclear and electronic 
eigenf unctions and the construction of the unperturbed 
eigenfunction and the basis functions ip n together 
with a discussion of the convergence properties is given 
by Taylor et al. §]. 

The electric dipole operator must also be transformed 
to the body-fixed axis. The necessary procedures are 
described by Lefebvre-Brion and Field [Q. For matrix 
elements of h ■ r connecting £+ states to £+ states 

(AT,A = 0|n-r|AT + l,A = 0) 

= [(Af + 1)/3] 1 / 2 (A = 0|2|A = 0) (15) 

and connecting £+ states to ±1„ states 

(N, A = 0|n-r|7V + l,A= ±1) = =f[(JV + 2)/3] 1/2 

x (A = 0\2- 1/2 {x =f iy)\A = ±1), (16) 

where r = (x,y,z). The calculation of (0|z|0) and 
(0\x^fiy\±l) in prolate spheroidal coordinates is straight- 
forward. 



Calculations of S(p) were carried out with basis sets 
comprised of electronic and vibrational functions P,[l0||. 
The converged values of S(Q) and S(— 1) obtained us- 
ing 121 electronic and 11 vibrational basis functions are 
given in Table [|. The convergence of the sum rules 
with basis set size is approximately logarithmic. Er- 
rors were determined for each sum S(p) by finding A 
and c such that Ae~ cn is the difference between the 
values obtained with basis sets of sizes n x n x n and 
(n + 1) x (n + 1) x [n + 1). The total error given in Ta- 
ble I for each entry is AJ^ =n e~ ct = Ae~ cn [l - e" ]" 1 . 

The values of the calculated sums S(0) and S(— 1) 
agree with the exact values jll],[l2| to better than 2 parts 
in 10 8 . Table | also lists the values of 5(-2), S(-3), 
and S(— 4). We anticipate no loss of accuracy in eval- 
uating S(— 2) since the summation Eq. (^) is stationary 
with respect to first order errors. The corresponding val- 
ues of the dipole polarizabilities ad are given in Table || 
and Fig. 0. The sums S(— 3) and S(— 4) are related to 
quantities occurring in the determination of the polar- 
izabilities [fl^3|-[L5|. S(— 3) enters in the combination 
B e = |5(— 3) — j£jC , where C is the scalar quadrupole 
polarizability. With C = 23.99 for H 2 + and 23.24 for 
D 2 + @, we predict that B 6 = 7.77 for H 2 + and 7.24 for 
D2 + . The empirical value for H2" 1 " derived by Jacobson 
et al. is 7.8(5). 

TABLE I. Nonrelativistic evaluation of the sum S(p), 
Eq. (ph, for H 2 + and D 2 +. 



p 


H 2 + 


D 2 + 





1.000 000 0(1) 


1.000 000 0(2) 


-1 


1.653 650 96(2) 


1.635 744 78(6) 


-2 


3.167 000 94(1) 


3.071 152 0(2) 


-3 


6.780 745 959(7) 


6.375 365 3(3) 


-4 


15.889 406 225(5) 


14.325 799 4(6) 



Comparison of theoretical nonadiabatic values 
dipole polarizability for the ground states of 
with experimental values. The results from 
have been multiplied by the factor (1 + e) 2 . 



TABLE II. 
of the electric 
H 2 + and of D 2 
Refs. Q and | 



H 2 



D 2 



Ref. 



3.168 0" 1 



RJTUOTS 
-0.0001 

3.168 2(4) 
3.168 5 
3.168 3 
3.168 725 6(1) 
3.1681(7) 



3.067 r 



lU.UtJlb 
-0.0020 

3.0714(4) 
3.071 87 
3.071 78 
3.071 988 7(2) 
3.0712(7) 



5 1, variational 
2 1, finite element 
4|, artificial channel 
1 1, variational 
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polarizability (atomic units) 

FIG. 1. Polarizabilities of Hb + and L>2 + in their ground 
states. For each of the two calculations from the present work 
the error bar is within the vertical line crossing through the 
data point. 

Table [n] and Fig. [l] contain a comparison of our calcu- 
lated values of ad with experiment and with the results of 
other theoretical methods. We leave aside calculations of 
the polarizability corresponding to an electric field along 
the body-fixed axis p7|JT^ ] . Moss Q employed a varia- 
tional method and an artificial channel method, with a 
classical description of the rotation. We are able to re- 
produce his results with our procedure if we take N = 
for the intermediate states with the consequent neglect 
of E — II coupling, the error introduced by ignoring rota- 
tional coupling being accordingly one in the fourth dec- 
imal place in the calculated polarizability. The calcu- 
lations of Bhatia and Drachman (|] and Shertzer and 
Greene 0] make no approximations other than in the 
numerical applications of their methods and yield values 
consistent to within the precision they claim with our 
results. 

We have determined the non-relativistic electric dipolc 
polarizabilities of the lowest rotational state of Hi2 + and 
D2 + to a precision, we believe, of one part in 10 8 . We 
expect that relativistic corrections will enter at the level 
of one part in 10 5 based on known corrections for the 
hydrogen atom @. Other effects arising from the finite 
size of the nucleus and nuclear spin will be still smaller. 
A new analysis of the experimental data (!]] incorporating 
our values of the sum rules may yield improved estimates 
of other properties that enter the interpretation. 
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